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ABSTRACT 


We present Gemini Planet Imager (GPI) observations of AU Microscopii, a young 
M dwarf with an edge-on, dusty debris disk. Integral field spectroscopy and broadband 
imaging polarimetry were obtained during the commissioning of GPI. In our broadband 
imaging polarimetry observations, we detect the disk only in total intensity and find 
asymmetries in the morphology of the disk between the southeast and northwest sides. 

The southeast side of the disk exhibits a bump at l" (10 AU projected separation) that 
is three times more vertically extended and three times fainter in peak surface brightness 
than the northwest side at similar separations. This part of the disk is also vertically 
offset by 69±30 mas to the northeast at l” when compared to the established disk mid¬ 
plane and consistent with prior ALMA and Hubble Space Telescope /STIS observations. 
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We see hints that the southeast bump might be a result of detecting a horizontal sliver 
feature above the main disk that could be the disk backside. Alternatively when in¬ 
cluding the morphology of the northwest side, where the disk mid-plane is offset in the 
opposite direction ~50 mas between O'M and 1".2, the asymmetries suggest a warp-like 
feature. Using our integral field spectroscopy data to search for planets, we are 50% 
complete for ~4 Mj up planets at 4 AU. We detect a source, resolved only along the disk 
plane, that could either be a candidate planetary mass companion or a compact clump 
in the disk. 

Subject headings: circumstellar matter — planet-disk interactions — stars: individual 
(AU Mic) — techniques: high angular resolution — instrumentation: adaptive optics 
— methods: data analysis 


Introduction 


Debris disks are comprised of rocky and icy bodies, left over from planet formation. The 
planetesimals in these disks collide, fragment, and produce micron-sized grains that scatter optical 
and near-infrared radiation from the star, resulting in a detectable nebulosity. Grains are influenced 
by the gravitational pull of larger bodies in the system as well as radiative forces and corpuscular 
drag due to the host star. These effects can induce morphology in debris disks such as midplane 


warps, as in j3 Pictoris (Mouillet et al. 1997), bow structures as seen in HD 61005 (Hines et al. 


2007), and narrow, offset rings such as Fomalhaut (Kalas et al. 2005). 


The debris disk around AU Microscopii (AU Mic) is interesting for multiple reasons. At a 
distance of 9.9 pc, we are able to probe closer to the star than most systems. As a member of the 
/3 Pictoris moving group (Barrado y Navascues et al.||1999), its young age of 23 ± 3 Myr (Mamajek 


& Bell 2014) makes it a prime candidate for directly detecting young planets still hot from their 


recent formation. It is also one of the very few M dwarfs with a debris disk imaged in scattered 


light. Since its ground-based discovery by Kalas et al. (2004), this edge-on disk has been studied 


in scattered light using both adaptive optics and the Hubble Space Telescope ( 

Liu 

2004. 

Krist et 

al. 2005 

Metchev et al. 2005 

Graham et al.| 2007 

Fitzgerald et al. 2007 

Schneider et al.|2014 

), as 

well as in millimeter thermal emission with interferometric arrays ( 

Wilner et al. 

2012; 

MacGregor 


et al. 2013). From modeling these observations, it is believed AU Mic harbors a “birth ring” of 
planetesimals at ~40 AU, inwards of which is empty and outwards of which is comprised of micron¬ 


sized dust particles placed on barely bound orbits (Strubbe &; Chiang 2006; Augereau &: Beust 


2006). Within the disk, there are also local enhancements of surface brightness and deviations in 
the disk’s vertical structure (e.g. Fitzgerald et al. 2007; Schneider et al. 2014). These features are 


likely non-axisymmetric structures located at or exterior to the birth ring and seen in projection 
(Graham et al. 2007). 


The Gemini Planet Imager (GPI) is a high contrast imaging instrument on Gemini South 
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with a high-order adaptive optics (AO) system, a coronagraph to suppress starlight, and a science 
instrument that features both integral field spectroscopy (spectral mode) and broadband imaging 
polarimetry (polarimetry mode) (Macintosh et al.|2014 Larkin et al.|2014 Perrin et al.|2015). GPI 


gives us the opportunity to probe closer in and deeper than previous optical or near-infrared obser¬ 
vations of AU Mic. We use GPI in both spectral and polarimetry modes to study the morphology 
of the debris disk and search for planets close in to the star. 


Observations and Data Reduction 


AU Mic was observed as part of commissioning the GPI instrument. On 2014 May 12, twenty- 
seven 60 s spectral mode frames were taken in the K1 spectral band (1.90-2.19 ^m) at an airmass < 
1.01. Images of the sky and instrument thermal background were not taken with the observations, 
but thermal background data were taken earlier that night and used in the data reduction. Forty- 
four 60 s polarimetry mode observations were taken on 2014 May 15 in H -band at an airmass < 
1.015. The GPI waveplate was rotated by 22?5 between frames. AU Mic transits almost directly 
overhead at Cerro Pachon, so a total of 154° of field rotation was obtained in each data set for 


Angular Differential Imaging (ADI, Marois et al. 2006). For the H -band data, four frames were 


discarded due to anomalies associated with commissioning tests. 

The raw data consist of individual spectra in spectral mode and pairs of spots for the two 
orthogonal polarizations dispersed by the Wollaston prism in polarimetry mode. Each spectrum or 
spot pair corresponds to a spatial resolution element, or spaxel, in a datacube created by the GPI 
Data Reduction Pipeline (DRP, Perrin et al. 2014). The spectral mode data were dark subtracted, 
corrected for shifts in the position of individual spectra on the detector due to instrument flexure 
(Wolff et al.| 2014) with an R-band argon arc lamp exposure taken immediately before the observing 
sequence, wavelength calibrated with K1 -band argon arc lamp data, fixed for bad pixels in the 2-D 
data, assembled into a spectral data cube, thermal background subtracted, fixed for bad pixels in 
the 3-D datacube, and corrected for distortion. 


The reduction of the polarimetry mode data follows the steps described in Perrin et al. (2015). 


To create polarimetry datacubes where the third dimension contains two orthogonal polarizations, 


the data were dark subtracted, flexure corrected using a cross-correlation routine (Draper et al 


2014), fixed for bad pixels in the 2-D data, assembled into a datacube using a model of the point 


spread function of the spot pair, fixed for bad pixels in the 3-D datacube, and corrected for dis¬ 
tortion. The polarimetry datacubes were then collapsed in the polarization dimension to obtain 
2-D total intensity images. For our total intensity reductions, we do not construct the full Stokes 
datacube, but one is made with the GPI DRP to search for a polarized signal. 

To combine and photometrically calibrate our data, we used the GPI DRP to measure the four 
fiducial diffraction or “satellite” spots, which are centered on the occulted star and imprinted with 


its attenuated spectrum (Wang et al. 2014}). In spectral mode, the satellite spot fluxes were used 
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to calibrate the data photometrically relative to the occulted star. For each datacube, the location 
of the occulted star at each wavelength was found using a least squares fit to all of the satellite 
spots’ positions and the magnitude of the atmospheric differential refraction. In polarimetry mode, 
the satellite spot flux calibration is an ongoing task, so no flux-calibrated polarimetry mode data is 
presented here. The central star in the polarimetry data was located using a radon-transform-based 


algorithm (Wang et al. 2014 Pueyo et al. 2015). The precision on the star center is 0.7 mas and 


0.9 mas for spectral and polarimetry mode respectively (Wang et al. 2014). 


To remove the stellar point spread function (PSF), we used pyKLIP (Wang et al. 2015), a 


Python implementation of the Karhunen-Loeve Image Projection (KLIP) algorithm (Soummer et 


al. 

2012; 

Pueyo et al. 

2015 

). To distinguish astrophysical sources from the stellar PSF, we used 


ADI for the polarimetry mode data and both ADI and spectral differential imaging (SDI, Marois 
et al.]|2000) for the spectral mode data. Multiple PSF-subtracted images were made depending on 
the purpose of the reduction. For all the subtractions, we divided the PSF into annuli, divided 
each annulus into azimuthal sectors, and ran KLIP on each sector. 


The main parameters we adjusted were the number of modes used from the Karhunen-Loeve 
(KL) transform and an exclusion criteria for reference PSFs, similar to Ng in Lafreniere et al. (2007), 
defined by the number of pixels a hypothetical astrophysical source would move azimuthally, radially 
or both due to ADI/SDI. For the H -band total intensity data, we used three KL modes and a ten- 
pixel exclusion for the diffuse left side of the disk presented in Figure Qa) and four KL modes and 
a six-pixel exclusion for the more compact right side presented in Figure Qa). We also searched for 
fainter features using ten KL modes and a three-pixel exclusion as shown in Figure [ljb) . For the 
K1 -band data, we used ten KL modes and a three pixel exclusion criteria to characterize a compact 
candidate source shown in Figure [l](c) and twenty KL modes and a 1.5 pixel exclusion criteria to 
filter out extended disk features and quantify our sensitivity to point sources in Section [4j 


Lafreniere et al. (2007 


The frames were rotated to orient the disk horizontally using a disk position angle (PA) of 
128?41 ( |MacG regor et al. 2013) derived from observations with Atacama Large Millimeter/submillimeter 
Array (ALMA) before being mean combined to produce the images in Figure [I] 

We did not detect any polarized signal (5% upper limit on the polarization fraction at 0'. / 5), 
consistent with a low expected polarization for forward scattered light at small projected angular 
separations (Graham et al. 2007) and reduced polarization sensitivity due to image persistence in 
the HAWAII 2-RG detector (Smith et al. 2008), for which calibration work in ongoing. 


3. Analysis of the Debris Disk 

All three reductions in Figure [T] show stark asymmetries between the southeast (SE) and 
northwest (NW) sides of the disk. To quantify these asymmetries, we used the conservative 77- 
band reduction in Figure [l](a) and fit Gaussian vertical profiles to the disk at various separations. 
We fit between (//4 and F/25. The inner radius is determined by the increasing noise at small angular 






































- 6 - 


separations and the outer radius is set by GPI’s finite field of view. Due to the relatively low signal- 
to-noise and ADI self-subtraction in the vertical wings of the disk, a Gaussian provides good fits 
with well-behaved residuals and performs as well as alternate profiles such as the Lorentzian used 


by Graham et al. (2007). The primary aim of the Gaussian fit was to measure the disk’s mid-plane 


vertical offset and full-width-half-maximum (FWHM). We also obtained the peak brightness of the 
disk as part of the fit and calculated the vertically integrated flux by summing over the flux within 
one FWHM centered about the measured disk mid-plane. 

Our reductions are subject to varying levels of disk self-subtraction, which biases our mea¬ 
surements. To correct for this effect, we repeated the analysis with model disks injected separately 
at one of five non-overlapping positions, resulting in ten sides for analysis. The disk model is a 
symmetric, edge-on, optically-thin ring parameterized by a power-law grain size distribution, a 
Henyey-Greenstein scattering function, and a Gaussian vertical profile. We adopted values from 


the porous water-ice model from Graham et al. (2007) and varied both the scale height and the 


surface brightness to accommodate the range of measured values. As GPI only sees the inner arc- 
second of the disk, the main contribution to the flux in the models is the front edge of the disk 
with scattering angles < 15°. The mean FWHM, peak brightness and integrated vertical flux we 
measured from the model disks were used to correct biases in our measured values as a function of 
separation. The scatter of all four quantities in the model disks were used as the uncertainties in 
our measurements. 


We present measurements of the mid-plane offset and FWHM as a function of separation 
with 1 -a error bars in Figure [2| The vertical mid-plane offset peaks around 1 " on the SE side 
with an offset from the average ALMA mid-plane of 69±30 mas to the northeast at 1" separation. 
The mid-plane offset is roughly constant between Cf/4 and T/2 on the NW side, with an offset of 
47±8 mas to the southwest at 1 ". The FWHM is roughly constant on the NW side and consistent 
with measurements at larger separations but triples on the SE side past (//8 to a peak FWHM 
of 570±56 mas. Due to unsubtracted residual PSF features, the disk brightness as a function of 
separation that we measure is not reliable. However, the overall brightness asymmetry between the 
two sides as seen in Figure [l] is real: on average, the peak surface brightness is three times larger 
on the NW side, but the vertically integrated brightnesses are consistent with being equal. 

From the aggressive reduction in Figure [ljb), we see a dark gap and a faint horizontal sliver 
on the SE side, above the main disk. This sliver is only at 3.1a significance when comparing the 
total flux of this feature to randomly summing up noise in the same manner. However, it emerges 
in successive reductions that progressively remove speckle noise and does not point radially toward 
the star, like residual speckles. We cannot rule out that the apparent morphology is influenced by 
KLIP self-subtraction. We also do not observe any indication of a corresponding feature on the 
NW side, which could be hidden underneath the over-subtraction of the brighter main belt of the 
disk. 








4. Sensitivity to Planets 


The achieved sensitivity for our K1 data was calculated following the procedure outlined by 
Mawet et al. (2014) and the 5cr-equivalent false positive threshold is plotted in Figure |3j Flux 
attenuation due to PSF subtraction was quantified by injecting and recovering the brightnesses of 
simulated planets. To translate our contrast curve to limits on possible planets in the system, we 


ran a Monte-Carlo analysis as described by Nielsen Sz Close (2010) to determine the completeness 


of our data. Since GPI is not sensitive to planets at large orbital distances due to its restricted held 
of view, we also included the limits from the Near-Infrared Coronagraphic Imager (NICI, Wahhaj 
et al. 2013) in Figure [3j With the improved inner working angle, the GPI data are 50% complete 
for 4 Mj up planets with semi-major axes of 4 AU assuming COND models (Baraffe et al.| 2003). 

We did find one compact K1 candidate, marked in Figure [l])c), at a separation of 544=L4 mas, 
a PA of 130?2=L0?9, and a K1 contrast of 1.14=0.3 xl0~ 5 . Although this source is detected with 5-a 
confidence in AT, there is no corresponding source at H. Because polarimetry mode flux calibration 
is still ongoing, we cannot currently quantify the constraints placed by our H-band data. Accord¬ 
ing to COND models (Baraffe et al. 2003), this source would correspond to a 2-3 Mj up planet. 
However, inspection of Figure [I|c) reveals that the clump is extended along the disk. The FWHM 
perpendicular to the disk is ~70 mas and consistent with being unresolved when compared to the 
satellite spot PSF; parallel to the disk the measured FWHM is ~160 mas giving an intrinsic size 
of ~140 mas (1.4 AU), suggestive of a local enhancement in the disk similar to those reported at 
larger projected separations (e.g., Liu 2004 Fitzgerald et al. 2007 Schneider et al. 2014). We 


are confident in the measured elongation because it is significantly larger than a residual speckle, 
which, given the K1 bandwidth, would subtend about 100 mas at this angular separation. 


Discussion 


The AU Mic debris disk has significant asymmetries in the morphology and surface brightness 
between the SE and NW sides of the disk. In Figure [4j we compare our observations with 1.3- 
nnn continuum data from ALMA (MacGregor et ah] 2013) and optical imaging using the using 


the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (Schneider et al. 


2014). The ALMA data probe larger planetesimals insensitive to the variable radiation pressure 


and stellar wind drag of AU Mic, making it easier to disentangle dynamical and radiative forces on 
the debris disk. STIS, like GPI, traces scattered light from small grains, but with a larger held of 


view. 


The disk mid-plane offsets, FWHM enhancement, and brightness asymmetries we find agree 


well with the published STIS results from Schneider et al. (2014). The most prominent common 
feature is the faint, extended, and vertically-offset bump on the SE side at 1". This feature is also 
continuous in the STIS data where it is seen to taper back down near 2" from the star. Since 
the SE bump is also consistent with the most prominent undulation in the ALMA contours, it is 










































plausible that the feature is associated with gravitational perturbations of the parent bodies. The 
statistical significance of this mm feature is questioned by MacGregor et al. (2013); however, the 
correspondence between ALMA and shorter wavelength observations lends credence to its reality. 
Assuming the mm grains are not displaced from the birth ring, the correlation between scattered 
light and mm-emission suggests that the SE bump is located in the ~40 AU birth ring and is just 
seen projected at ~10 AU separation. As both the STIS and ALMA data were taken within four 
years of the GPI observations, we cannot yet distinguish between orbital motions at 10 and 40 AU. 


The compact K1 candidate in our data could be another clump of disk material or a planet 
with a circumplanetary disk. The candidate must be a transient feature unless the object is self- 
gravitating. If the projected size of this source is representative of its radial extent (da), Keplerian 
shear at an orbital separation (a) of 40 AU will completely erase this feature in only (2/3 ){a/5a) = 
19(a/40 AU)(1.4 AU/da) orbits. However, a 1.4 AU Hill sphere at 40 AU is established by a 
0.08(M*/0.6M©) Mj up object, which is well below our detection limits. Follow up observations are 
required though for verification and to determine common proper motion. 


Ignoring the possible detection of the sliver in Figure [IJb), the vertical mid-plane offset varia¬ 
tions on the two sides of the disk in opposite directions hint at a vertical warp in the disk. Applying 
the analytical model from Dawson et al. (2011) to the AU Mic system, any inclined planet between 
an Earth and Neptune mass could induce the warp depending on the orbital separation and dynam¬ 
ical timescale. The candidate companion we see at a projected separation of 5.5 AU can be slightly 
lower than an Earth mass and still be responsible for the warp. The FWHM and flux asymmetries 
in the disk could be due to dynamical heating from planets or parent bodies in the ring, which can 
puff up the disk vertically but keep the total amount of starlight scattered constant. 


If the sliver feature is real, it may be the backside of the disk (the bright main belt would be 
the front side). In this case, the variations in the disk could be explained by an inclined ring with 
azimuthally varying brightness due to clumps in the disk, like those seen at further separations 
(Fitzgerald et al. 2007). Essentially, on the SE side, the backside is visible due to a local brightness 
enhancement that is absent on the NW side. However, without a corresponding detection on the 
northwest, the backside interpretation is uncertain. 


Our GPI observations show that there remain many unanswered questions regarding this exten¬ 
sively studied debris disk. The stark asymmetries between the two sides of the disk and correlation 
with millimeter observations hint at gravitational interactions due to perturbers. With deeper 
imaging by new near-infrared high contrast imaging instruments and higher resolution data in the 
millimeter by ALMA, we may soon have a much better understanding of the physics in the system. 
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Fig. 1.— Reduced images of the AU Mic disk. The images are rotated so that the disk is horizon¬ 
tal. Images (a) & (b) are from the 77-band total intensity broadband observations. Image (a) is 
optimized for the bright main belt of the disk. The azimuthal structure outside of the horizontal 
disk plane is due to unsubtracted speckles. Image (b) is used to look for fainter features. The curly 
brace indicates the position of the sliver feature we detect. Image (c) is from the Kl -band spectral 
mode observations. The arrow indicates the candidate compact source that we detect. 
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Fig. 2.— Measurements of the morphology of the disk in H -band with the negative and positive 
separations for the SE and NW sides respectively. At each separation, we plot the disk mid¬ 
plane vertical offset and FWHM with 1-cr error bars as blue squares. The faded circles represent 
the FWHM before correcting for biases. Adjacent points are correlated so that the plot is well 
sampled. 
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Fig. 3.— (Top) Sensitivity to point sources as a function of projected separation using the Kl -band 
spectral mode data. (Bottom) Completeness for planets as a function of semi-major axis and mass 
in the AU Mic system based on contrast curves from NICI (Wahhaj et al. 2013) and GPI (this 


paper). Completeness is calculated using a Monte Carlo method described by Nielsen & Close 


(2010) where planets with random orbits are generated, the contrast curves determine whether 
they are detected, and COND models (Baraffe et al. 2003) are used to convert from luminosity to 


mass. 
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Fig. 4.— Comparison of GPI data with previous observations. GPI H -band data from Figure [lja) 


is shown in the center with contours from 1.3-mm ALMA observations ( 

MacGregor et al. 2013| 

overlaid on top using a linear scaling and an image from STIS ( 

Schneider et al. 

2014 

) in the 


background. The southeast (left) bump between —1.5 ,/ and — 1" in the disk is apparent in all three 


observations. The mid-plane offset is also continuous between the STIS and GPI data. 








